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(54) Precede de modelisation d'une surface et dispositif pour sa mise en GBuvre. 

fe7) L'invention concerne un precede d'obtention d'une 
moaelisation d'une surface. 

II est caracterise par les etapes consistant a obtenir des 
mesures de donnees geometriques relatives a des points 
specifiques de la surface, effectuer un maillage de la sur- 
face passant par lesdits points, memoriser a une adresse 
specifique a chaque nceud du maillage les coordonnees du 
nceud, le nombre de satellites du noeud, une information 
d'acces aux adresses desdits satellites et par suite aux in- 
formations y relatives et des donnees geometriques even- 
tuellement associees audit noeud; pour chaque noeud, defi- 
nir un indice de rugosite local obtenu a parti r d'une somme 
ponderee des coordonnees actuelles du noeud et de ses 
satellites, definir la somme d'un indice de rugosite global, 
constitue par ('addition de tous les indices de rugosite lo- 
caux, et d'un indice global de violation desdites donnees 
geometriques, ajuster de facon iterative les coordonn6es 
des noeuds indetermines, en utilisant a chaque fois I'addi- 
tion d'une combinaison pond6ree des coordonnees actuel- 
les de voisins du noeud et d'une combinaison des donnees 
geometriques associees au noeud, pour amener ladite 
somme a un minimum, et creer une representation de la 
surface a partir des coordonnees ajust6es. 
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La presente invention concerne d'une fagon 
g6n6rale 1' investigation dans des corps ^ trois 
dimensions, et concerne en particulier un nouveau precede 
pour l'obtention d'une representation d'une surface 
situ6e au sein de ce corps k partir d'un ensemble limits 
de donn6es g6ometriques connues relatives k ladite 
surface. 

L 1 investigation dans des corps tridimensionnels 
1Q est une preoccupation majeure tout particuliferement dans 
les domaines de la geologie, de la g6ophysique et de la 
biologie. 

Ainsi, par exemple dans le domaine de la 
geophysique, on souhaite obtenir des representations 
aussi fiddles que possibles de surfaces situees par 

15 exemple k 1' interface entre deux zones de natures ou de 
propri6t6s diff6rentes, k partir de donnfees obtenues lors 
de campagnes de prospection, voire en cours 
d' exploitation d'une richesse souterraine. 

En effet, l'efficacite de 1 1 investigation du corps 

20 tridimensionnel, et en particulier de la recherche de 
p6trole ou autre richesse souterraine, est conditionn6e 
par la fid61it6 avec laquelle on peut reconstituer et 
repr6senter ce genre de surface. 

Par ailleurs, dans le domaine medical et 

25 biologique, l'on sait actuellement obtenir, par des 
proced6s divers, des representations de coupes d'un corps 
vivant ou analogue, et l'on recherche 6galement k 
reconstituer et k representee de fa<?on 
tridimensionnelle, une surface situ6e k 1" interface entre 

30 deux milieux, et notamment les contours d'un organe ou 
autre. 

II existe k 1'heure actuelle un certain nombre de 
techniques de modeiisation permettant d' obtenir la 
representation d'une surface dans un corps k trois 
35 dimensions que l'on cherche k exploiter ou k traiter. On 
citera en particulier les interpolations de B6zier ou les 
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fonctions splines (on se r6f6rera pour plus de details 
par exemple k l'ouvrage Geometric Modelling, M. E. 
Mortenson, Ed. John Wiley, 1985) . 

Cependant, ces techniques connues sont toutes tr6s 

5 mal adapt6es k la prise en compte de donn6es h6t6rog6nes, 
k savoir des donn6es relatives k la fois aux coordonn6es 
de points de la surface, k des orientations impos6es de 
la surface, k des liens g6om6triques entre deux points 
distants, etc. En outre, ces techniques sont inadapt6es k 

10 une trop grande dispersion des surfaces ou encore k des 
anomalies de surfaces telles que des plis, des cassures, 
ou encore des surfaces discontinues. Plus pr6cis6ment, 
dans certains cas, les fonctions mises en oeuvre dans ces 
techniques soit ne convergent pas, soit n'admettent 

15 aucune solution, soit admettent au contraire une 
plurality de solutions. 

Enfin les techniques connues n* off rent pas la 
possibility de prendre en compte la notion de plus ou 
moins grande certitude avec laquelle certaines donn6es de 

20 la surface sont connues. 

La pr6sente invention vise k pallier les 
inconv6nients des proc6d6s de la technique ant6rieure et 
k proposer un proc6d6 qui puisse prendre en compte des 
informations g6om§triques non homog6nes et/ou tr6s 

25 disperses, ainsi que les anomalies pr6cit§es et 
d'autres, en off rant k chaque fois une mod61isation 
unique et sans Equivoque de la surface. 

Un autre objet de 1' invention est de proposer un 
proc6d6 qui puisse s'appuyer sur des donn<§es g6om6triques 

30 en prenant en compte simultan6ment un degr6 de certitude 
ou de precision li6 k ces donn^es. 

A cet effet, la pr§sente invention concerne tout 
d'abord un proc§d6 d'obtention d'une mod61isation d'une 
surface representative notamment de 1 ' interface entre 

35 deux zones de natures ou de propri6t6s diff6rentes dans 
un corps k trois dimensions tel qu'une. formation 
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g§ologique ou un corps vivant, caract6ris6 en ce qu'il 
comprend les 6tapes consistant k: 

- obtenir au moyen d'appareils de mesure un 
ensemble de donn§es g6om6triques relatives k la surface 

5 et associ£es respectivement k des points donn6s de ladite 
surface; 

- effectuer un maillage de la surface de telle 
sorte que 1' ensemble desdits points donn6s soit un sous- 
ensemble des noeuds du maillage; 

10 - m6moriser k une adresse donn6e d'une m6moire, 

sp6cifique k chaque noeud du maillage f les informations 
suivantes: 

* les coordonnees du noeud consid6r6, 

* le nombre de noeuds satellites du noeud 
15 consid6r§, 

* une information permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives & ces derniers, 

* le cas 6ch6ant des donn6es g6om£triques 
20 associ£es audit noeud consid£r6 r 

- pour chaque noeud du maillage, d6finir un indice 
de rugosity local obtenu a partir d'une somme pond6ree 
des coordonnees actuelles du noeud et de ses satellites, 

- d6finir la somme d'un indice de rugosite global, 
25 constitu6 par 1* addition des indices de rugosite locaux 

associ6s k chaque noeud, et d'un indice global de 
violation desdites donn6es g6om6triques ; 

- ajuster de fa<?on iterative les coordonnees de 
chacun des noeuds dont on ne connalt pas les coordonnees 

30 precises, en utilisant & chaque fois I 1 addition d'une 
combinaison pond6r6e des coordonnees actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donn6es g6om6triques associ6es 
audit noeud, de manidre k amener ladite somme k un 

oe minimum; 
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- creer une representation de la surface & partir 
des coordonn6es ajust6es de chacun des noeuds. 

L * invention concerne 6galement un dispositif pour 
l'obtention d'une modelisation d'une surface 
representative notamment de 1' interface entre deux zones 
de natures ou de propri6t6s diff6rentes dans un corps k 
trois dimensions tel qu'une formation g£ologique ou un 
corps d'§tre vivant f caract6ris6 en ce qu'il comprend : 

- des moyens de mesure un ensemble de donn£es 
g6om6triques relatives & la surface et associ6es 
respectivement & des points donnas de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donn6s soit un sous-ensemble des noeuds du maillage; 

- une m6moire contenant & des adresses donn6es, 
specif iques a chaque noeud du maillage, les informations 
suivantes : 

* les coordonn6es du noeud consid6r6, 

* le nombre de noeuds satellites du noeud 
consid6r6 f 

* une information permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donn6es g6ora6triques 
associ£es audit noeud consid6r£, 

- des moyens de calcul pour ajuster de fagon 
iterative les coordonn6es de chacun des noeuds dont on ne 
connalt pas les coordonn£es pr6cises, en utilisant A 
chaque fois 1 1 addition d'une combinaison ponder 6e des 
coordonn6es actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
donn£es g6om6triques associ§es audit noeud, de mani£re & 
amener a un minimum une somme d'un indice de rugosite 
global, constitu§ par 1' addition d f indices de rugosite 
locaux respectivement associ6s h chacun des noeuds et 
chacun obtenu & partir d'une somme pond6r6e des 



2652180 



coordonnees actuelles dudit noeud et de ses satellites, 
et d'un indice global de violation desdites donnees 
g§om6triques ; et 

- des moyens pour cr6er une representation de la 
5 surface k partir des coordonnees ajustees de chacun des 
noeuds • 

D'autres aspects, buts et avantages de la pr6sente 
invention apparaitront mieux k la lecture de la 
description d6taill6e suivante d'une forme de realisation 
10 pr6f6r6e de celle-ci, donn6e k titre d' exemple non 
limitatif et f aite en reference aux dessins annexes , sur 
lesquels : 

la figure 1 est un exemple de graphe incorporant 
un maillage utilise pour modeiiser une surface, ce graphe 

15 etant defini par 1* ensemble de tous les cotes de facettes 
dont les sommets constituent 1' ensemble Q d6crit plus 
loin, et illustre egalement des filches symbolisant des 
d6placements de noeuds du maillage effectu£s 
interactivement par 1 ' utilisateur pour fagonner la 

2 Q surface; 

les figure 2(a) k 2(d) illustrent quatre cas de 
controle de la forme d'une surface k l'aide de quatre 
types differents de contraintes vectorielles ; les 



figures 2(a) k 2(c) et connus par leur orthogonality A un 
vecteur V donn6 sur la figure 2(d); 

la figure 3 illustre un exemple de surface 
geologique complexe (d6me de sel) modeiisee k l'aide du 
proc6d6 de la pr6sente invention; 

la figure 4 illustre un exemple de surface 
biologique complexe (cerveau d'embryon) modeiis6e k 
1 ' aide du proc6d£ de la presente invention k partir de 
donn6es ponctuelles obtenues par une succession de 
coupes; 




sont supposes connus directement sur les 



35 



la fxgure 5 illustre 1 'organisation de la memoire 
associ6e k un atome A(k) de noyau ; 
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la figure 6 illustre un premier exemple de 
memorisation d'une orbite d'un atome donne, dans un 
tableau de nk cases; 

la figure 7 illustre un second exemple de 
5 memorisation d'une orbite d'un atome donne, dans nu 
tableaux cha£n6s de deux cases chacun; 

la figure 8 illustre la memorisation d'une donn6e 
g6om6trique ou contrainte; 

la figure 9 illustre la maniere de m6moriser un 

10 

ensemble de donn6es g6om6triques ou contraintes relatives 
& un atome A(k) donne ; 

la figure 10 illustre la memorisation dans un 
tableau de 1' ensemble des atomes relatifs h une surface 
S; et 

15 la figure 11 illustre 1 ' enchalnement de 1 1 acc£s 

audit ensemble d ! atomes. 

On indiquera tout d'abord que, d'une figure & 
1' autre, des elements ou parties identiques ou similaires 
sont d6sign6s par les m§mes num£ros de reference. 

20 En outre, on va indiquer ci-dessous un certain 

nombre de publications dont les contenus respectifs sont 
a consid6rer comme inclus dans la pr6sente description. 

(1) A method of Bivariate interpolation and smooth 
surface fitting based on local procedures, H. AKIMA , J. 

25 ACM 17.1, 1974; 

(2) Shape reconstruction from planar cross-sec- 
tions, J»D. BOISSONNAT, ACM Transactions on Graphics, 
vol. No. 2, 1986; 

(3) Machine contouring using minimum curvature, 
30 I.C. BRIGGS, Geophysics 17,1, 1974; 

(4) Triangular Bernstein-B£zier patches, G. FARIN, 
Computer-aided geometric design, 1986; 

(5) Primitives for the manipulation of general 
subdivisions and the computation of Voronoi diagrams, L. 

35 GUIBAS et J. STOLFI, ACM Transactions on Gr aphis , vol. 
No. 2, 1985; 
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(6) Three-dimensional graphic display of 
disconnected bodies, J-L. MALLET, Mathematical Geology, 
vol. 20, No. 8, 1988; 

(7) Geometrical modelling and geostatistics , J.L. 
5 MALLET, Proceedings of the Third International 

Geostatistics Congress, Kluwer Academic Publishers, 1989? 

(8) Discrete Smooth Interpolation, J.L. MALLET, 
ACM Transactions on Graphics, April issue, 1989; 

(9) Geometric Modelling, M.E. MORTENSON, John 

!0 Wiley, 1989. 

Par ailleurs, la description qui suit d'une forme 
de realisation prSferee de "l' invention est elle-m§me 
suivie d'une annexe exposant une methode de 
reconstitution d'une surface d<§riv6e de la m6thode 
expos6e dans le document (8) et sur laquelle est bas6 le 

15 

proc6d6 de la prfesente invention. 

1. DESCRIPTION DU PROCEDE DE CODAGE ET D 1 AJUSTEMENT 
DE SURFACES COMPLEXES SELON L 1 INVENTION 



1.1 Introduction 

20 1.1.1 Objet de la description 

On va decrire ci-dessous un procede industriel de codage et 
d'ajustement de surfaces complexes rencontrees par exempie lors de la moderation 
d 'interfaces entre zones de natures ou de proprieties differentes de couches geologiques 
ou de corps d'organismes biologiques. 
25 Le resultat de ce codage peut etre transcrit sous forme d'images et/ou de maquettes 
en materiel solide (par exempie plastique) a l'aide d'un ordinateur graphique ou de 
tout autre procede. En geologie, par exempie, ces images et/ou maquettes sont ensuite 
utilisees pour optimiser I'exploration, Sexploitation et la gestion des ressources du sous- 
sol telles que : 

30 • petrole 

• minerai 

• eau 

• etc ... 

Le procede decrit peut etre decompose en deux sous-procedes compiementaires : 
35 1. un sous-procede de codage d'une surface 
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2. un sous-procede d'ajustement d'une surface aux donnees disponibles. En geologie, 
par exemple, ces donnees peuvent etre : 

• la position (exacte ou approchee) de points sur une coupe de la surface par 
un plan 

5 • la position (exacte ou approchee) de ['intersection d'un forage avec la surface 

geologique correspondant a la separation entre deux couches 

• le plan tangent (exact ou approche) en un point donne de la surface corre- 
spondant a la separation entre deux couches 

• le rejet (exact ou approche) d'une faille (=cassure) affectant la surface cor- 
1 0 respondant a la separation entre deux couches 

• etc ... 

Le sous-procede de codage est con^u de fa$on a optimiser les performances du procede 
d'ajustement base sur une variante de la methode DSI (cf ref. [8]) specialement adaptee 
a la moderation geometrique des surfaces. Cette variante de la methode DSI utilisee 
1 5 ici n'a encore jamais ete publiee et est presentee en annexe. 

1.1.2 Decomposition d'une surface S en facettes 

Soit 5 une surface a modeliser. Comme le montre la Figure 1 et ainsi qu'il a ete 
suggere dans les references [6] et [7], celle-ci sera supposee approchee par un ensemble 
de facettes triangulaires et/ou polygonales. Ces facettes sont caracterisees par : 

20 • La connaissance des coordonnees de leurs sommets constituant un ensemble fini de 
N points . . . , tp N } regulierement repartis sur 5. Ces sommets sont numerotes 
de 1 jusqu'a iV et les coordonnees du fc eme sommet tp k seront notees (y?f, <p y ki <p z k ) 

• La donnee de leurs cotes constitues par des liens entre couples de sommets (£,-, (p^) 

L'ensemble des cotes des triangles et/ou polygones composant une surface S constitue 
25 un graphe G dont les noeuds sont precisement les sommets des triangles et/ou poly- 
gones. Pour cette raison, et afin d'eviter des periphrases inutiles, nous utiliserons de 
fa^on equivaiente le mot "noeud" a la place de "sommet de triangle et/ou polygone" : 

noeud = { sommet de triangle et/ou polygone } 

1.1.3 Definitions preliminaires 
30 Ensemble des satellites A(ife) 

Soit <p k le fc^ me noeud constituant une surface 5. Nous appellerons ensemble des 
satellites de $ k et nous noterons A(k) Tensemble des noeuds differents de (p k tels 
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<l ue soit un cote de triangle ou de polygone. En d'autres termes, 1 ensemble 

A(fc) est tel que : 



10 



1) = % ^} 

merae 
gone 



ff k sont les extremites d'un 
cote de triangle ou de poly- 



Comme indique dans la relation ci-dessus f nous conviendrons de noter n k le nombre de 
satellites rattaches au noeud (p k . 

Cette notion d'orbite introduce dans les articles [6] et [7] n'a jusqu'a present pas 
fait 1'objet de procede de codage particulier; une partie iirpartante de l 1 invention 
consists a proposer un procede de codage pour A(Jt). 



Voisinage N(k) 

Soit tp k le * eme noeud d'une surface S. Nous appellerons voisinage de $ k et nous 
noterons N{k) l'ensemble des noeuds constitues par la reunion de (p k et de ses satellites. 
En d'autres termes, on a : 

15 = U A(Jb) 

Nous appellerons voisin du noeud $ k tout noeud appartenant a l'ensembie JV(Jfc). En 
d'autres termes, on a : 

<? y e N(k) <j=> {(p j = Voisin de <p k ] 

Cette notion de voisinage introduite dans les articles [7] et [8] n'a jusqu'a present pas 
20 fait 1'objet de procede de codage particulier; une partie iitportante de 1' invention 
consiste a proposer un procede de codage pour N(fc). 

Atome A(k) et noyau associe 

Soit <p k le * eme noeud d'une surface 5. Nous appellerons atome de noyau (p k et nous 
noterons A(k) le couple (<p k , A*(*)) constitue par la juxtaposition (et non la reunion) : 

25 • du noeud <p k 

• de l'ensembie A*(t) des adresses ou des codes . . . , A*(;„J} permettant 

de retrouver dans une memoire l'ensembie des atomes {A(h), .... A(j nk )} corre- 
spondants dont les noyaux sont les satellites de (p k 
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En d'autres termes, on a : 

n k — nombre de satellites de (p k 

avec < A*(j) ~ Adresse ou code d'acces de l'atome A(j) 

correspondant au noyau du i eme satellite 
de <p k 

Cette notion d'atome introduite dans Tarticle [6] n'a jusqu'a present pas fait I'objet de 
procede de codage particulier; une partie importance de l 1 invention consiste 
a proposer un procede de codage pour A(k), 

Contraintes et types de contraintes 

10 Nous appellerons contrainte relative a un atome A(Jb), toute condition relative a la 
position du noeud ip k correspondant au noyau de cet atome. Par exemple, en geologie on 
peut citer les contraintes suivantes qui sont d'une importance capitale pour Tajustement 
d'une surface aux donnees sures ou approchees observees ou mesurees par des procedes 
manuels ou geophysiques : 

• Les coordonnees {<P k ,<p k i<p%) du noeud (p k sont connues de fa^on exacte. Une 
telle contrainte est dite de type noeud de controle et est utilisee en geologie 
par exemple pour coder Tintersection de la surface modelisee avec un forage. 

• Les coordonnees (sp£»{p£i?£) du noeud $ k sont connues de fa^on approchee avec 
un coefficient de certitude donne. Une telle contrainte est dite de type noeud 
de controle flou et est utilisee pour coder la position approximative d'un noeud 
de la surface modelisee supposee mesuree de fa$on approximative par un procede 
donne. Par exemple, 

- en geologie, les donnees sismiques sont de ce type, 

- en biologie, les donnees de coupes d'organes observees sous microscope ou 
25 par echographie ou encore par scanner sont de ce type. 

• Le vecteur (<p x —<p k ) joignant le noeud $ k a un autre noeud <p x est connu de faqon 
approchee avec un coefficient de certitude donne. Une telle contrainte est dite 
de type lien vectoriel flou et est utilisee en geologie par exemple pour coder le 
rejet d'une faille, le noeud £ fc etant sur un bord de la faille et le noeud (p x etant 

3Q alors sur le bord oppose (voir Figure 2 ). 

• Le vecteur V orthogonal au plan tangent a la surface S au point correspondant 
au noeud tp k est connu de fagon approchee avec un coefficient de certitude donne. 



15 



20 
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Une telle contrainte est dite de type normale floue et est utilisee en geologie par 
exemple pour coder le pendage des couches mesure par des procedes manuels ou 
geophysiques (pendagemetrie). 

• etc ... 

5 Les coefficients de certitude mentionnes ci-dessus sont des nombres positifs dont la 
valeur est censee representer la confiance que l'on peut avoir a a priori" dans les donnees 
qui lui sont attachees. La prise en compte des contraintes floues presentees ci-dessus 
constitue Tune des caracteristiques.de l'invention^st pr^sentfe en detail 
dans I'annexe (voir annexe - section 2.5 et 2.7.6). 

10 Convention 

Afin de simplifier les notations, on convient d'identifier les numeros des noeuds d'une 
surface S avec les noeuds eux-memes. II s'ensuit que Ton a : 

• A(fc) = ensemble des satellites {Vh*'^*Vj nk } de tp k ou ensemble des numeros 
{h i • • • . Jnjk } correspondants 

15 • N(k) = ensemble des voisins {$ k , $ h , . . . , <p jnk } de <p k ou ensemble des numeros 
{k t ji , . . . , jn k } correspondents 

1.1.4 Etat de l'art 
Les methodes classiques 

On trouvera dans, la reference [7] une description tres complete de Tetat de Tart concer- 
20 nant la modelisation de surfaces complexes. En particulier, dans cet article, on expose 
les raisons pour lesquelles les methodes classiques utilisees en cartographie automatique 
et en Conception Assistee par Ordinateur (CAO) ne sont pas adaptees a la modelisation 
de telles surfaces comme on aurait pu le croire a priori. A la fin de cet article [7], on 
introduit brievement les fondements de la methode DSI (quelquefois appelee DSI/DSA) 
25 permettant de surmonter les difficultes inherantes aux methodes classiques. 

La methode DSI 

Cette methode est basee sur la notion de "rugosite" d'une surface S composee de 
facettes triangulares et/ou polygonales. Tout d'abord, la notion de rugosite R{jp\k) au 
noeud <p k d'une surface S est definie a Taide de la formule suivante (Voir ref. [7], [8] et 
30 Annexe) : 

Les coefficients positifs ou nuls {v a (k)} intervenant dans la formule ci-dessus sont des 
ponderateurs choisis par Tutilisateur charge de modeliser la surface. Parmi Tinfinite de 
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choix possibles, l'un des plus interessants appelle "ponderation harmonique" consiste a 
poser : 

V c(u\-S si <* = k 

v W-\ 1 si a€A(Jfe) 

Les rugosites locales ainsi definies sont ensuite combinees en une rugosite globale R(jp) 
definie comme suit : 

s(£) = £K*)-W) 
k 

Les coefficients positifs ou nuls jx(fc) intervenant dans la formule ci-dessus servent a 
moduler la contribution des rugosites locales R(<p\k) a la rugosite- globale Si Ton 

decide par exemple d'utiliser une ponderation uniforme, il suffit de poser : 

/i(*) = l VJt 

Parmi les autres choix possibles pour les coefficients on peut decider d'utiliser par 
exemple les ponderateurs suivants ou m est une constante donnee : 




1 si <p k est un point de donnee 
m > 1 si (p k n'est pas un point de donnee 



Le principe de la methode DSI consiste a ajuster la position de chacun des noeuds 
<p k de la surface S modelisee de telle sorte que celle-ci soit la plus lisse possible tout 
en respectant les donnees et contraintes regissant sa forme, L'article [8J mentionne en 
bibliographie est consacre a une presentation mathematique approfondie de la methode 
DSI et propose une technique de minimisation de la rugosite globale respectant les 
donnees et contraintes. 



Version interactive de la methode DSI 

La methode de minimisation de la rugosite globale de la methode DSI proposee dans 
Tarticle [8] est mal adaptee a une utilisation interactive sur calcuiateur electronique 
(ordinateur ou microprocesseur specialise), c'est pourquoi nous proposons en annexe 
une autre technique de minimisation originate sur laqueile est basee l 1 invention. 
Ainsi qu'il a et' dit precedemment, Tobjet de l'inventionest precisement de proposer 
un procede de codage et un procede d'ajustement de la surface S suivant cette nouvelle 
version de la methode DSI. 
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1.2 Exemples damages d'objets 

codes et ajustes par notre procede 

1.2.1 Exemple d'image de surface geologique 
obtenue par notre procede 

5 La recherche du petrole est generaiement effectuee de faqon indirecte en ce sens que 
plutot que de rechercher le petrole, on cherche des couches geologiques reservoir sus- 
ceptibles de le contenir ou de le capturer; par exemple, la deformation des couches 
de sel donne naissance a ce que Ton appelle des "domes de sel" qui sont consideres 
comme generant d'exceilents pieges a petrole. Pour pouvoir localiser avec precision la 

10 position du petrole il est necessaire de connaitre la forme exacte des couches "pieges" 
et malheureusement dans de nombreux cas il n'existait jusqu'a ce jour aucun procede 
pratique permettant de les modeliser de fagon satisfaisante; le cas des domes de sel est 
a ce titre considere comme Tun des plus complexes et pour prouver Tefficacite de notre 
procede de codage et ajustement nous presentons sur la Figure 3 Pimage d'une telle 

15 surface obtenue avec celui-ci. 

1.2.2 Exemple d'image de surface biologique 
obtenue par notre procede 

Les procedes d'acquisition de donnees tels que scanner ou echographie utilises en 
medecine sont matheraatiquement identiques aux procedes d'acquisition de donnees 

20 sismiques utilises en prospection petroliere. Par exemple, Techographie est identique 
a la sismique reflexion et le fonctionnement des scanners est base sur la tomographic 
tout comme la tomographic sismique. Des lors, il n'est pas etonnant que notre procede 
de codage et ajustement des surfaces soit egalement utilisable pour modeliser la peau 
d'organes bioiogiques et nous proposons sur la Figure 4 le modele du cerveau d'un 

25 embryon humain de 7 millimetres de long obtenu a Paide de notre procede de codage et 
d' ajustement. II est anoter qu r un aussi petit organe (le oerveau fait 1 milli- 
metre de long) ne pouvait jusqu f a present etre observe que sous forme de 
coupes microscqpiques. 



1.3 Procede de codage d'une surface S 

30 1.3.1 Notion de memoire et d'adresse de memoire 

Nous appellerons "memoire" tout dispositif de calculatrice electronique qui enregistre 
des donnees, des codes ou resultats partiels qu'U restitue au moment opportun pour les 
faire intervenir dans la suite des operations. De plus, nous appellerons "adresse" d'une 
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memoire, tout code permettant de retrouver Templacement physique de la memoire a 
Tinterieur de i'ordinateur. 

1.3.2 Codage d'un atome A(k) 

Nous proposons de coder et stocker chaque atome A(Jfc) dans un ensemble de memoires 
5 a i'aide du procede suivant ou Tune de ses variantes : 

Proced6 de codage 

Comme le suggere la Figure 5 , dans une zone memoire contigue commen$ant a l'adresse 
A*(k), on stocke de faqon consecutive les informations suivantes relatives a l'atome A(k) 
correspondant au noeud <p k ; 

10 1. noyau = 

• soit les coordonnees ({Pfi^i^t) du noeud $ k 

• soit l'adresse ou le code ou toute autre information permettant de retrouver 
les coordonnees (^f,^,^) du noeud $ k 

2. nb-sat = 

15 • soit le nombre njt de satellites lies du noeud $ k 

• soit l'adresse ou le code ou toute autre information permettant de retrouver 
le nombre de satellites lies du noeud $ k 

3. sat = 

Adresse d'une zone memoire contenant les adresses ou les codes d'acces 
20 . . . , A*(j nk )} permettant de retrouver les n k atomes {A(ji) t . . . , A(; n J} 

dont les noyaux {tp h , , . . , $ jn ^ } sont les satellites du noeud (p k 

4. const = 

Codage de Tensemble des contraintes (cf section 1.1.3) pouvant affecter le noeud 
<p h . Ce codage sera presente en details dans la section 1.3.3. 

25 5. info = 

Zone memoire de taille variable contenant zero, une ou plusieurs informations 
complementaires concernant Patome A(k). Ces informations sont specifiques a 
^application particuliere envisagee; par exempie en prospection petroliere, si A(k) 
est un atome correspondant a la surface separant deux couches d et C 2 , alors il 
30 peut etre interessant de stocker dans le champ info des informations telles que : 

• liste des proprietes physiques (porosites, vitesse des ondes sismiques, ...) de 
la couche Ci au noeud $ k correspondant a Patome A{k) 
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• liste des proprieties physiques (porosites, vitesse des ondes sismiques, ...) de 
la couche C2 au noeud (p k correspondant a l'atome A(k) 

• liste des proprieties geologiques (facies geologique, presence de petrole, ...) 
de la couche C\ au noeud <p k correspondant a Tatome A{k) 

5 • liste des proprietes geologiques (facies geologique, presence de petrole, ...) 

de la couche C2 au noeud <p k correspondant a Fatome A{k) 

• etc .. 
Variante 1 

La zone memoire dans laquelle sont stockees les adresses ou les codes d'acces 
10 {A*(ji) t . . . , A*(jn k )} permettant de retrouver les n k atomes {A(;'i), . . . , A(j nfe )} dont 
les noyaux sont les satellites {p^ , . . . , <p^ } du noeud ip k peut etre structuree de deux 
facpns : 

1. La premiere fac,on illustree par la Figure 6 consiste a utiliser un tableau de n k 
memoires consecutives contenant les adresses {A*(ji) t . . . , A*(j nk )} 

15 2. La seconde fac,on illustree par la Figure 7 consiste a utiliser n k tableaux cora- 
portant chacun 2 memoires consecutives telles que pour le or eme de ces tableaux : 

• la premiere memoire contienne Padresse A*(j a ) de I'atome correspondant au 
a* m * sateUite de I'atome A(k) 

• la seconde memoire contienne 1'adresse du (a + l)* me tableau. Si a est egal 
20 au nombre n* de satellites (nb-sat=rifc), alors cette seconde memoire est 

soit inutilisee, soit contient un code indiquant que Ton a atteint le dernier 
satellite de A(k) 

Variante 2 

II est possible de ne pas stocker le champ nb-sat mais dans ce cas, si la liste des 
25 adresses ou codes d'acces {A*(,?i), . . . , A*(j nk )} est codee sous forme d'un seul tableau 
(voir variante 1), il faut rajouter un emplacement memoire A*(n* + 1) dans lequel on 
place un code indiquant que la liste A*(Jfc) est terminee. 

Variante 3 

Le champ sat peut etre choisi suffisamment grand de facpn a y stocker directement les 
30 adresses ou codes d'acces {A*(ji), . „ , A*(j nk )} 
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Variante 4 



Un grand nombre de variantes equivalentes sont possibles suivant la taille du champ 
info et de la facon dont on' redecoupe celui-ci. II est toutefois bien entendu que ces 
variantes n'ont aucune importance pour ['utilisation de l'atome A{k) en connection avec 
l'algorithme d'ajustement base sur la methode DSI decrite en annexe. 



Variante 5 



On peut rajouter d'autres champs a l'atome A{k) decrit ci-dessus mais ceux-ci n'ont 
aucune importance pour l'utilisation de l'atome A(Jfc) en connection avec l'algorithme 
d'ajustement base sur la methode DSI decrit en annexe. 

Variante 6 

Quelles que soient les variantes, la seule chose importante pour ^implementation efficace 
de Palgorithme DSI est la suivante : 

Ayant l'adresse ou le code d'acces A*[k) d'un atome A(Jfc), U est necessaire 
d'acceder de la fa$on la plus directe possible aux atomes dont les noyaux 
sont les satellites du noyau de A{k). 

L'object essentiel du pcocede de codage de A(k) dans une memoire decrit ci-dessus ainsi 
que toutes ses variantes est precisement de permettre Faeces direct aux satellites. 

1.3.3 Codage des contraintes 

Chaque contrainte relative a un atome A(k) sera stockee dans une zone memoire 
specifique dont la description est proposee dans ce qui suit : 

Proc4d6 de stockage d'une contrainte 

Le codage d'une contrainte particuliere relative a un atome donne est effectue dans une 
zone memoire comprenant les deux champs contrainte et next comrae 1'indique la 
Figure 8 . Ces deux champs sont utilises de la fagon suivante : 

• contrainte = 
Zone memoire contenant : 

- soit les informations relatives* la contrainte codee (voir section plus bas) 

- soit l'adresse ou le code d'acces d'une zone memoire contenant effectivement 
les informations relatives a la contrainte codee 

Parmi les informations relatives a la contrainte codee, doit toujours figurer un 
code permettant de reconnaitre le type de la contrainte codee (cf section 1.1.3). 
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• next = 

zone memoire contenant : 

- soit un code indiquant que la contrainte consideree est la derniere des con- 
traintes affectant Patome A(k). 

5 - soit Padresse d'une zone memoire contenant la contrainte suivante (voir sec- 

tion suivante). 

Procede de stockage de toutes les contraintes affectant un atome 

V 

Le codage d'un atome A(k) decrit dans la section 1.3.2 contient un champ const que 
nous utiliserons pour stocker : 

10 • soit un code indiquant que Patorae A(k) n'est affecte par aucune contrainte 

• soit Padresse d'une zone memoire contenant la premiere contrainte affectant 
Patome A(k) 

Si Padresse de la n eme contrainte relative a Patome A(k) est stockee dans le champ next 
de la (n- l) eme contrainte (n > 1), alors, comme le suggere la Figure 9 , on est assure 
1 5 que ie procede de codage decrit permet d'attacher a Patome A(fc) toutes les contraintes 
qui lui sont relatives. Ce mode de codage permet de rajouter ou de supprimer facilement 
des contraintes sans remettre en question Porganisation des autres donnees. 

Exemples d' informations relatives k une contrainte donn£e 

Dans la section 1.1.3 nous avons presente a titre d'exemple trois types de contraintes 
20 tres importants pour le codage et Pajustement des surfaces complexes rencontrees en 
biologie ou geologie. A chacun de ces types de contraintes correspondent des informa- 
tions qui doivent etre associees au champ contrainte decrit au paragraphe precedent 
et schematise sur la Figure 9 : 

• Pour une contrainte de type noeud de controle, en plus du type, ii est necessaire 
25 de stocker au minimum les donnees suivantes : 

- les coordonnees ((fi t (p y ks <p z k ) approchees du noeud <p k correspondant au noyau 
de Patome A{k) auquel est attachee la contrainte 

- le coefficient de certitude (nombre positif) mesurant la confiance que Pon 
peut avoir dans les coordonnees (y^,^,^). 

30 • Pour une contrainte de type lien vectoriel, en plus du type, il est necessaire de 
stocker au minimum les donnees suivantes : 

- Padresse du noeud <p x auquel est lie le noeud ip k correspondant au noyau de 
Patome A{k) auquel est attachee la contrainte 
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- les trois composantes approchees du vecteur {jp x - <p k ) joignant le noeud $ k 
a au noeud (p x 

- le coefficient de certitude (nombre positif) mesurant la confiance que Ton 
peut avoir dans les composantes approximatives du vecteur (£ A - (p k ) . 

5 • Pour une contrainte de type normale, en plus du type, il est necessaire de stocker 

au minimum les donnees suivantes : 

- les trois composantes approchees du vecteur V orthogonal a la surface S au 
noeud <p k correspondant au noyau de Tatome A{k) auquei est attachee la 
contrainte 

10 - le coefficient de certitude (nombre positif) mesurant la confiance que Ton 

peut avoir dans les composantes approximatives du vecteur V . 

La presentation detaillee de ces trois types de contraintes est proposee dans les sections 
2.5 et 2.7.6 de i'annexe. 

Variante 

15 Bien que Ton n'ait parle dans ce qui precede que d'un seul coefficient de certitude 
pour les composantes en x, y et z de chaque contrainte, on pourra affecter a ces 3 
composantes 3 coefficients de certitude differents si necessaire. 

1.3.4 Codage d'une surface S 

Afin d'optimiser la programmation de la variante de la methode DSI presentee en 
20 annexe, nous proposons de coder une surface 5 composee de facettes triangulares 
et/ou polygonales sous forme d'un ensemble de N atomes dont les noyaux (noeuds) 
sont les sommets de ces triangles et/ou polygones. Les donnees et codes correspondant 
a chacun de ces atomes {A(l), . . . , A(N)} peuvent etre stockes : 

• soit dans un tableau comme le suggere la Figure 10 

25 • soit dans le champ atom de zones memoires chainees entre elles par leur champ 
next comme le suggere la Figure 11. Dans ce cas, le champ next associe a 
Tatome numero k doit contenir Tadresse ou le code d'acces de i'atome suivant 
numerote (k + 1); le champ next associe au dernier atome A(N) doit contenir un 
code indiquant que le dernier atome a ete atteint. 

30 La deuxieme solution est la plus interessante si Ton souhaite pouvoir ajouter ou sup- 
primer facilement des atomes, par contre elle consomme davantage de place memoire 
que la premiere. 
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1.3.5 Variantes 

Pout des commodites de presentation, les zones memoires utilisees pour coder des 
informations ont ete decoupees en champs, et a chacun de ces champs a ete attribue : 

• un nom propre 

5 • une place dans la zone memoire correspondante 

A est bien evident que Ie procede de codage decrit dans ce brevet ne depend pas des noms 
ni de l'ordre de ces champs a I'interieur des zones memoire utilisees pour le codage. Par 
aiUeurs. il est toujours possible de rajouter des champs a des fins specifiques a certaines 
applications sans que cela modifie le procede de stockage. 



10 1.4 Ajustement d'une surface S 
1.4.1 Introduction 

Etant donne une surface S codee suivant le procede decrit dans la section 1.3, nous 
nous proposons dans ce qui suit de decrire un procede base sur ce codage et permettant 
d'ajuster de fagon optimaie les coordonnees de certains des noeuds .... $ N \ et ceci 
1 5 de sorte que la surface S : 

• soit la plus lisse possible au sens du critere de rugosite globale de la methode DSI 
(voir references [7], [8] et annexe ainsi que la section 1.1.4) 

• respecte le mieux possible les informations et contraintes connues avec plus ou 
raoins de certitude concernant la forme de cette surface ( voir sections 1.1.1 et 

20 1-1.3) 

Pour ce faire, on utilisera une variante originate de la methode DSI presentee en annexe 
specialement concue pour utiliser le codage decrit dans la section 1.3. Compte tenu des 
informations disponibles observees ou mesurees, cette methode ajuste successivement 
et de facon iterative les trois coordonnees de chaque noeud $ a non fixe 

25 suivant une formule du type : 

nouvelle valeur de ^ = /£( anciennes valeurs de <p% , . . . , <fF N ) 
nouvelle valeur de ^ = f%( anciennes valeurs de <p{, . . . , <p y N ) 
nouvelle valeur de ^ = /*( ; anciennes valeurs de ^f, <p%) 

Plus precisement, si l'on se reporte a l'annexe presentant la variante utUisee de la 
30 methode DSI, ces fonctions d'ajustement /£(), et £() sont deduites de la forme 
locale des equations DSI au point numero a et ont la forme suivante (voir sections 2.4.1, 
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2.4.2 et 2.7.3 de l'annexe) 



10 



/«(•••) = 



-]k- { E {/<(*)*"(*)• E A*)-v£} 

~Wy { £ «*>*■<*>■ E A*)-*g> 

t ' 

-j?.- { E M*) *•(*)■ E A*)-**) 

i * 



0±a 



0*a 



15 



avec : < 



= E M(i)(v a (fe)) 2 

fc€AT(a) 



La nature et le role des coefficients de ponderation {v a (k)} et ont ete indiques 
20 section 1.1.4 et sont detailles dans I'article [8] ainsi qu'en annexe. Les valeurs des 
termes {If , IT. If }, {if. if*, if*} et des termes {Q^^TAD dependent des 
types et des valeurs des contraintes affectant le noeud ip a et leur formulation exacte est 
presentee en annexe. 

Dans le cas particuUer ou les coefficients de ponderation {v a (k)} correspondent 
25 a la ponderation harmonique presentee section 1.1.4 et en annexe, les fonctions de 
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ponderation f£Q, f*Q et /*() se simplified et prennent la forme suivante : 
«-) =]^-[ /<«) • |A(«)| • C £ rf) 

0 V fc€A(a) 7 

- E K*)-{(E -4 

- E{rr a -<3r> ] 

• |A(a)| ■ ( £ tf) 

E k*)-{(e ^)-iawiv4 
Ear -on ] 

K") • |A(a)| • ( £ tf) 

W(a) ' 

E m*h(e ^-iawi-^} 
Er-«n ] 



avec : < 



M a =M«)-|A(a)| 2 + £ K*) 

*SA(or) 

1 5 Nous proposons dans ce qui suit un precede de calcul de ces fonctions base sur le codage 
presente au paragraphe 1.3: 

1.4.2 Notations 

Afin de simplifier la description des precedes de calculs presences dans ce qui suit, il 
est commode d'introduire les notations et definitions suivantes : 



2652180 



22 



• Par definition, nous appellerons bloc cooperations, toute liste d'operations 
precedee d'une accolade ouvrante et terminee par une accolade fermante u }" . 
Les operations figurant dans un bloc d'operations peuvent etre des operations 
elementaires et/ou d'autres blocs d'operations. 

5 • Par definition, nous appellerons memoire de travail toute memoire utilisee pour 
stocker des resultats intermediates au cours d'un procede de calcul. 

• Par definition, nous appellerons "'chargement" de 1 'expression e dans une memoire 
m l'operation consistant a evaluer la valeur de e puis a ranger ie resultat dans la 
memoire m. Par convention, une telle operation sera notee de la fa^on suivante : 

10 m <— e 

• Pour toute memoire m utilisee pour stocker une valeur numerique, nous ap- 
pellerons "incrementation" de m par l'increment i l'operation consistant a calculer 
la valeur m+ i puis a ranger le resultat dans la memoire m. Par convention, une 
telle operation sera notee de la fagon suivante : 

15 m «— m + i 

Dans l'operation ci-dessus, 1'increment i peut etre : 

- soit une constante 

- soit le contenu d'une memoire 

- soit une expression arithmetique 

20 Dans l'operation ci-dessus, l'increment i peut etre : 

• soit une constante 

• soit le contenu d'une memoire 



• soit une expression arithmetique 
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1.4.3 Procede d'ajustement dans le cas de ponderateurs harmoniques 



Dans le cas ou les coefficients {v a (k)} sont les ponderateurs harmoniques presentes dans 
la section 1.1.4 et dans Tannexe, les valeurs /£, /g et des fonctions d'ajustement 
/50i faO et faO relatives au noeud <p a peuvent etre calculees de facon efficace a l'aide 
5 du procede suivant base sur le codage de S : 



10 



0) AUouer les memoires de travail suivantes 



15 



f r f y f* 

J at J J a 



20 



1) EfFectuer les initialisations suivantes (rordre n'a aucune importance) 



25 



<j% = a * = a % = 0 



30 



2) Determiner le nombre n a de satellites de A(a). Ce nombre peut etre obtenu 
a partir du champ nb-sat de A(a) 

3) Soit A(a) i'ensemble des satellites de l'atome A(a) auxquels on accede par 
le champ sat de cet atome. Pour chaque satellite k € A(a), repeter les operations 
du bloc suivant : 

{ 
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15 



20 



3.1) Determiner le nombre de satellites de A(k). Ce nombrje 
peut etre obtenu a partir du champ nb jsat de A(k) 

3.2) Determiner les coordonnees du noeud (p k de 
1'atome A(k). Ces coordonnees peuvent etre obtenues a partir 
du champ noeud de A{k). 

3.3) A partir des informations contenues dans A{k) % calculer ou 
deter miner la valeur du ponderateur fi(k) (voir section 1.1.4 et 
annexe). 

3.4) Effectuer les operations suivantes (i'ordre n'a aucune impor- 
tance). 

*n 4 -« r n+* tf (*) 



s£ = si = 0 



3.5) Soit A(fc) Tensemble des satellites de 1'atome A(k) auquel 
on accede par le champ sat de cet atome. Pour chaque satellite 
0 € A(fc), 0 £ a effectuer les operations du bloc suivant : 

{ 

3.5.2) Determiner les coordonnees actuelles (<p£ip£ s pj) du 
noyau fip de i'atome A(0). Ces coordonnees peuvent etre 
obtenues a partir du champ noeud de A(f3). 

3.5.3) Effectuer les 3 operations suivantes (I'ordre n'a au- 
cune importance) : 



25 



30 



} 



3.6) Effectuer les 3 operations suivantes (Pordre n'a aucune im- 
portance) : 



} 



ft 
fa 
fa 



-£ + /<(*)x(s»-n t xtf) 



35 



4) Effectuer les operations suivantes : 



/<(a) x n a x n a + 



5) Soit C(a) l'ensemble des contraintes attachees a l'atome A(a) auquel on 
accede par le champ const de cet atome. Pour chaque contrainte c € C(a) repeter 



25 
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les operations du bloc suivant : 
{ 

5.1) Determiner le type de la contrainte c 

5.2) En fonction du type de la contrainte c, effectuer les deux 
operations suivantes : 

5.2.1) Extraire les informations specifiques de la contrainte c 

5.2.2) En fonction des informations extraites en (5.2.1), cai- 
culer les valeurs des expressions suivantes (voir sections 2.5 et 
2.7 de l'annexe) : 

pro pyor -r\zct 
4 i » k i » A i 

or. or, or 

5.3) Effectuer les operations suivantes (I'ordre n'a aucune impor- 
tance) 

M*' -Af** +7 f« 

n - n +(rr- or «) 



} 



6) Effectuer les operations suivantes (I'ordre n'a aucune importance) : 
fl - (M*) xn a xv* v - fi)/M? 

25 A la fin de ce procede de calcul, les memoires de travail /£, fg et /* contiennent les 
valeurs des fonctions d'ajustement /=(), /£() et /*() utilisees pour ajuster les valeurs 
des coordonnees du uoeud <p a : 

30 *4 /« 

Le procede d'ajustement decrit ci-dessus ne doit etre applique qu'aux noeuds dont la 
position est non fixee, c'est-a-dire aux noeuds qui ne sont pas des noeuds de controle 
(voir section 1.1.3). 



1.4.4 Variantes 
35 Variante 1 

La mise en oeuvre du procede d'ajustement (base sur la methode DSI) presentee au 
paragraphe precedent dans le cadre de ponderateurs harmoniques peut facilement etre 
adaptee dans le cas de ponderateurs quelconques. La structure globale du procede reste 
la meme, seule change 1 'incrementation des differ entes memoires de travail allouees a 
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l'etape (0); la fagon d'incrementer ces memoires de travail depend des ponderateurs 
choisis. 

Par exemple, si nous choisissons les ponderateurs {v a (k)} suivants, 



= k 

: € A(*) 



10 



15 



a . -1 si a-. 

alors on verifie facilement que les fonctions d'ajustement associees sont : 
£(-) =7^-f /•(«)•( E rf)/|A(«)| 

- E • {( E ^)/iA(i)i 2 -<P%/\m\\ 

0** 

- £{rf° - op} ] 

/*(*)•( E rf)/|A(ar)| 

V Jt€A(a) 7 

E m(*)-{(E ^)MWI 2 -^/|A(*)l} 

Jt€A(a) U /J€A(*) 7 J 

2{ir-QP> 1 

/*(«)•( E ^)/|A(a)| 

V *6A(a) 7 

E *»(*)•{( E ^)/|A(fc)[ 2 -^/|A(fc)|l 

*€A(a) '"Vac*) 7 J 



avec 



f Mr = iw„+E7f a 

= + Elf* 
M> = /*(<*)+ E W)M(*)I 2 } 

fc6A(o<) 

20 Le procede de calcul des valeurs f%, f% et /* prises par ces fonctions au noeud (p a est 
alors quasiment identique a celui decrit pour les ponderateurs harmoniques : 
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0) AUouer les memoires de travail suivantes : 



/<*> fat fa 



10 



15 



1) Effectuer les initialisations suivantes (I'ordre n'a aucune importance) 



a% = <r% = <r£ = 0 



20 



2) Determiner le nombre n a de satellites de A(a). Ce nombre peut etre obtenu 
a partir du champ nb-sat de A(a) 

25 3) Soit A(a) Tensemble des satellites de Patome A(a) auxquels on accede par 
le champ sat de cet atome. Pour chaque satellite k € A(a), repeter les operations 
du bloc suivant : 

{ 



28 
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3.1) Determiner le nombre n h de satellites de A(k). Ce nombre 
peat etre obtenu a partir du champ nb«sat de A(k). Effectuer 
ensuite l'operation suivante : 

n\ «- n k x n k 

3.2) Determiner les coordonnees {<p ki <p k ,<p k ) du noeud <p k de 
I'atome A(k). Ces coordonnees peuvent etre obtenues a partir 
du champ noeud de A(k). 

3.3) A partir des informations contenues dans A(k), calculer ou 
determiner la valeur du ponderateur /x(jfc) (voir section 1.1.4 et 
annexe). 

3.4) Effectuer les operations suivantes (Pordre n'a aucune impor- 
tance). 

5^=4=5* =0 

3.5) Soit A(fc) l'ensembie des satellites de I'atome A(k) auquel 
on accede par ie champ sat de cet atome. Pour chaque satellite 
0 6 A(£), P £cc effectuer les operations du bloc suivant : 

{ 

3.5.2) Determiner les coordonnees actuelles (^ 9 ^,^) du 
noyau <p p de I'atome A(ft). Ces coordonnees peuvent etre 
obtenues a partir du champ noeud de A{0) . 

3.5.3) Effectuer les 3 operations suivantes (I'ordre n'a au- 
cune importance) : 

3.6) Effectuer les 3 operations suivantes (Pordre n'a aucune im- 
portance) : 

& - fS +/*(*) * {s%/n\ - 9 %/n h ) 
fZ^f y a + M)x(synl-<pl/n k ) 
/^/5 + Kt)x(4/nl-^/n fc ) 
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4) Effectuer les operations suivantes : 

M*' «- M a 



10 



15 



20 



25 



5) Soit C(a) l'ensemble des coatraintes attachees a 1'atome A(a) auquel on 
accede par le champ const de cet atome. Pour chaque contrainte c G C(a) repeter 
les operations du bloc suivant : 

{ 

5.1) Determiner le type de la contrainte c 

5.2) En fonction du type de la contrainte c, effectuer les deux 
operations suivantes : 

5.2.1) Extraire les informations specifiques de la contrainte c 

5.2.2) En fonction des informations extraites en (5.2.1), cal- 
culer les valeurs des expressions suivantes (voir sections 2.5 et 
2.7 de i'annexe) : 

pra pya p za 
4 1 i i i L i 



5.3) Effectuer les operations suivantes (I'ordre n'a aucune impor- 
tance) : 

M?+yf a 

n-fi+ccr-on 



> 



m? 

Ml' 
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6) Effectuer les operations suivantes (I'ordre n'a aucune importance) : 

SI - 0*(a) x <r$/n a - fZy/M? 
fZ^W«)*<ryn a -fZ)/Mp> 

Variante 2 

Les pararaetres numeriques stockes dans le champ info des atomes peuvent egalement 
etre interpoles par la methode DSI (voir section 2.1.2 et 2.2 de 1'annexe). Dans ce 
cas, la structure du procede d'interpolation est quasiment identique a celle du procede 
d'ajustement de la surface presente ci-dessus; le seul changement consiste a ne conserver 
qu'une des trois fonctions d'ajustement normalement dediees aux trois coordonnees des 
noeuds et a calculer celle-ci en fonction du parametre interpole. 
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La pr6sente invention trouve application tout 
particuliSrement dans les domaines suivants : 

en g6ologie et en g6ophysique, elle permet de 
cr6er un module de surface, par exemple de construire une 
maquette solide en trois dimensions, sur laquelle on peut 
effectuer une simulation g6ophysique ou une simulation de 
1* exploitation de richesses souterraines : elle permet 
egalement d'obtenir un ensemble de donn<§es propres h 
initialiser un simulateur num6rique ou analogique de 
ph6nom6nes geophysiques ou d' exploitation de richesses; 
en outre, elle permet, & l'aide d'une station de travail 
graphique, d'afficher sur un <§cran ou d'imprimer sur un 
support des images permettant de controler les 
simulations ci-dessus; 

en biologie ou en m6decine, 1' invention permet de 
construire des modules en trois dimensions d'organes 
et/ou de proth^ses , ou encore de simuler les effets 
depurations en chirurgie plastique; de m§me, des images 

en deux dimensions affich6es ou imprim6es permettent ici 
encore de controler les operations ci-dessus. 

Bien entendu, la pr6sente invention n'est 
nullement limitee & la forme de realisation decrite ci- 
dessus et representee sur les dessins, mais l'homme de 
l'art saura y apporter toute variante ou modification 
conforme & son esprit. 
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ANNEXE 



2. Nouvelle version de la m£thode DSI adapt€e h 

la mod^lisation interactive des surfaces complexes. 



Dans Ies domaines de la Geologie et de la Biologie, il est frequent d'avoir a modeliser des 
surfaces complexes correspondant par exernple a des interfaces entre zones de nature 
ou de proprietes differentes. Les techniques classiques de modelisation basees sur Ies 
interpolations de Bezier ou sur les fonctions Splines (voir ref. [9]) sont tres mal adaptees 
a la prise en compte de ces donnees heterogknes. Nous proposons ici une approche 
differente basee sur la methode DSI (Discrete Smooth Interpolation - voir ref. [8]). Dans 
cette approche, les surfaces sont modelisees avec des facettes triangulares irregulieres 
dont Ies sommets doivent etre determines aGn de tenir compte du plus grand nombre 
de donnees heterogenes possible. 



2.1 Introduction : 

Presentation des problemes etudies 

2.1.1 Surface S definitions associees 

Comme le montre la Figure 1 , soit S une surface composee de facettes triangulares ad- 
jacentes plongee dans Tespace euclidien 3D (0\x, y t z). Cette surface n'est pas forcement 
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connexe ni fermee. Dans ce qui suit, nous utiliserons les notations suivantes : 

G = G{S) = ensemble des cotes des triangles 
i ft = ft(«S) = ensemble des sommets des triangles 
N = |ft| = nombre de points de ft 

G est un graphe dont les noeuds sont les points de ft et nous identifiers cet ensemble 
avec IWmble des indices {1, - r N} utilises pour les numeroter. Le voisinage iVTfc) 
d'un noeud fc € ft est defini par : 

N{k) = sous ensemble de ft constitue par les noeuds de ft qui 
peuvent etre atteints a partir du noeud fc en s(fc) pas 
1 0 au plus ou s(fc) est un nombre positif donne. 

Dans ce qui suit, nous supposerons que ces voisinages iV(fc) possedent la propriete de 
symetrie suivante : 

Les voisinages ainsi definis engendrent une topoiogie sur G qui peut etre consideree 
1 5 comme une approximation de la topoiogie de S elle meme. 

2.1.2 Probleme 1 : Estimation d'une fonction <p deflnie sur ft 

Supposons que la position de chaque noeud fc € ft est connue et soit tp(k) une fonction 
definie pour chaque noeud € ft et connue seulemeht sur un sous-ensemble de ft : 



L = ensemble des noeuds 1 6 ft oil <p(l) = <p t est connu 
I = ensemble des noeuds i e ft ou tp(i) = tp. est inconnu 
= ft-I 



Generalement, L est different de ft et il y a une infinite de fonctions <p{k) definies sur ft 
et interpolant les valeurs {<p(l) = ^ : f 6 I}. Notre objectif est de selectionner parmi 
toutes ces fonctions celle qui minimise un critere donne mesurant : 

25 • la rugosite globale de la fonction 

• l'ecart entre la fonction tp(-) et un certain nombre de contraintes et donnees qu'elle 
doit satisfaire 

2.1.3 Probleme 2 : Ajustement de la surface S 

Dans ce cas, seuls les points d'un s6us ensemble L de noeuds € ft ont une position 
30 connue £(fc) = fi k : 

L = ensemble de noeuds / g ft ou <p(l) = ft est connu 
/ = ensemble des noeuds igft oii = ft est inconnu 

= ft-L 
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Comme dans le cas du premier probleme, L est generalement different de f2 et il y a 
une infinite de surfaces S interpolant les points {$(1) = (p ( : I E L}. Notre objectif est 
de selectionner parmi toutes ces surfaces, celle qui minimise un critere donne 
mesurant : 

5 • la rugosite globale de la surface S 

• l'ecart entre la surface S et un certain nombre de contraintes et donnees qu'eile 
doit satisfaire 

Compte tenu de l'independance des composantes (^(Jfc), ^ y (fe), <f z {k)) de chaque vecteur 
£(fc), il est facile d'utiliser la solution du premier probleme pour resoudre le second en 
1 o posant simpiement : 

arm = + 

Dans le cas tres particulier ou Tensemble Q correspond aux noeuds d'une grille rectan- 
guiaire posee sur la surface, on peut songer a utiiiser une methode d'interpolation du 
type spline, Briggs ou Bezier (Voir ref. [9], [1] et [3]); malheureusement ces methodes 
1 5 sont inappliquables dans le cas des grilles irregulieres et ne permettent pas de tenir 
compte des contraintes presentees aux paragraphes 2.2.3, 2.5 et 2.7.2. 

2,2 Interpolation 

<Tune fonction (p(k) definie sur 

Dans cette section, nous introduisons brievement la methode DSI (Discrete Smooth 
20 Interpolation) presentee dans l'article [8]. 

2.2.1 Notations 

Dans ce qui suit, nous noterons <p une matrice coionne de taille N telle que : 

25 Notre probleme ne depend pas de la methode de numerotation des noeuds de la grille. 
Afin de simplifier les notations, nous supposerons que les elements de la matrice tp 
ont ete eventuellement permutes de telle fagon que <p puisse etre partagee en deux 
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sous-matrices <p t et <p L telles que : 



<Pl 



<Pr = 



<Pl = 



tin J 
Vim J 



avec i a € I V a 



avec IffEL V /? 



De plus, nous noterons || • \\ D une semi-norme associee a une matrice carree semi-definie 
positive [D] de taille (N x N) de telle facpn que, pour toute matrice colonne X de taille 
10 iV on ait : 

\\X\\l = X'-[D]-X 

2.2.2 Definition d'un critere de rugosite globale R((p) 

Soit R(<p\k) le critere de rugosite locale defini par la formule suivante ou les {v a (k)} 
sont des coefficients de ponder ation positife, negatifs ou nuls donn£s : 

a€iV(Jt) 

R((p\k) peut etre utilise pour construire un critere de rugosite globale Rfa) defini de la 
fagon suivante oii est une fonction de ponderation donnee definie sur ft : 

*(*0 = £m(*) 
k 

R((p) est completement defini par la donnee des coefficients {v a (k)} et fi{k) et il est 
20 facile de verifier (voir ref. [8]) qu'ii existe toujours une matrice (N x N) symetrique 
semi-definie positive [W] telle que : 

W) = /.[P7(fe)].^ 

2.2.3 Prise en compte de contraintes iineaires au sens des moindres 
carres 

25 Considerons une matrice carree [AJ de taille (AT x N) et une matrice colonne B { de 
taille AT definissant ia contrainte lineaire suivante : 

Etant donnee une matrice semi-definie positive [D{] } le symbole figurant dans la 
formule ci-dessus a la signification suivante : 



30 



I [Ai] -<P- B { Wj). aussi petit que possible 
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S'il existe plusieurs conditions de ce type, nous proposons de mesurer le degre de 
violation de ces conttaintes par le critere p(<p) defini comme suit : 

M = Ell W "P-*i \\ 2 D{ 



10 



15 



20 



2.2.4 Solution du probleme 

Parmi toutes les fonctions <p(k) definies sur et interpolant les donnees {y>(/) = p ( • 
, nous en choisiroas une qui minimise le critere suivant ; 

R?(<p) = R(<p) + p( v ) 
En developpant R*(<p), on obtient : 

R*(<p) = ^-[W*]-? - 2-<p t -[Q) + C 
'[W] = £ [AiY • [£>,] • [A { ] + [W] 

i 

avec \ Q =E Wf'&A-Bf 
C = E II B > f Di 

i 

La partition de la matrice <p induit une partition similaire des matrices [W*] et Q 
utilisees pour definir R*(<p) : 



^ = 



V>L 



[W*] 



- \ Wi 



LI W* L L 



; Q = 



Qi 
Ql 



(DSI) 



La condition dR*(<p)/d<p r = [0] impUque 1'equation suivante appelee "equation DSF 
et caracterisant toutes les fonctions {<p{k) : k 6 Q) solution de notre probleme : 

.avec:^ = -\Wf L ] ■ <p L + Q , 

2.3 Unicite de la solution de Pequation DSI 

Definition 

L'ensemble L des noeuds ou <p est connu sera dit consistant relativement a R*(<p) si 
25 chaque composante connexe du graphe G contient au moins un noeud appartenant a 
L. 
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theoreme 



Si le critere de rugosite globale R(<p) est tel que : 



5 



1) L est consistant relativeraent a R(*p) 

2) p( k ) >0 VJfc 

3) v a {k) >0 V {a £ fc, a € N(k)} 

4) «*(*) < - £ v°{k)£Q 



a6W(fc) 



aiors, i'equation DSI basee sur R{<p) admet une solution unique. 
Demonstration 

1 0 Le theoreme ci-dessus ne presente qu'un interet theorique et n'a aucune incidence sur 
le precede decrit dans le present brevet, e'est pourquoi nous ne donnerons pas ici sa 
demonstration. 

2.4 Forme locale de Pequation DSI 

Dans ce qui suit, piutot que de resoudre directement I'equation DSI, nous proposons 
15 une approche iterative originale permettant d'eviter le calcul et le stockage de [W*]. 
Cette approche iterative est precisement celle qui est utilisee dans le present brevet 

Afin de simplifier les notations, nous supposerons dans ce qui suit que les matrices 
semi-definies positives intervenant dans I'equation DSI sont diagonales. De plus, 
lorsqu'il y a des contraintes lineaires du type [A { ] • <p ~ B it nous utiliserons les 
20 notations suivantes : 




2.4.1 Calcul de dB?(<p)/d<p a 



Compte tenu de la definition de iZ(<p|fc), on peut ecrire Timplication suivante : 



= E E 



Pin 
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dR(<p\k) 
d<P a 



' 2 • (v Q (fc)) 2 • <p a + 2 • • £ t^(fe) • ^ si a G .V(jfc) 



* 0 sinon 



On deduit de cette relation que 



5 1 dRj<p) _ 1^ dR(<p\k) _ 1 ^ dR{<p\k) 

10 = E /»(*) ("°W) J + E E 



Compte tenu de la structure diagonale de [£),], il est facile de verifier que le a* me 
element de d || [A { ] <p-Bi \\ 2 D . /df a est tel que : 



15 



20 = *»«-E a** (4 a ) 2 + z{m*d?< • E 4^} - q? 

k k ^ 0J:a ' 

Compte tenu de la definition de R*(<p), on en deduit que : 



+ E M*)v a (*)- E A*)-^} 
+ EEf4^-E^} 

30 " ZQ? 
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La solution Vt est telle que dR*(<p)/d<fi a = 0, done la a^ e composante <p a de <p doit 
satisfaire l'equation DSI(a) suivante : 

= -jj*- { E E **(*)• l*> 



DSI(a) 



avec : < 



M °= E A*) (« a (fc)) 2 
if =Ea-*(4 q ) 2 

i? = E{^ fe -E^} 

Of = E A; a Di k Bi 
k 

Nous proposons d'appeler cette equation /o me /oco/e rfe /'egaofeon DSI zm noeud a. 
2.4.2 Proposition pour un algorithme iteratif 

15 La forme locale de l'equation DSI presentee dans la section precedente suggere un 
algorithme evident pour estimer la solution v . Par exemple, a Iteration (n + 1), la 
a composante y4 n+ > de la solution <pW doit satisfaire l'equation DSIfe) ce qui 
permet de proposer le procede iteratif suivant : 



20 



25 



soifc / l'ensemble des indices de noeuds ou (p a est inconnu 
soit tf> une solution approchee initiale 
tant_que( il est necessaire de faire de nouvelles iterations ) 
{ 

pour_tout( agl) 
{ 

M$ = {■■■} 

} *--*•{•••> 

} 



I PeUt 16 VOiri CCt al S orithn * simple n'utilise pas explicitement la 

matnee [W*] mtervenant dans l'equation DSI mais il requiert de recalculer un grand 
nombre de fois des produifcs du type {««(*) . qui sont pt ecisement ceux utilises 

30 pour construire [W]. 
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En fait, si la solution approchee initiate est proche de ia solution exacte, alors il 
n'est pas necessaire d'effectuer beaucoup d'iterations et le surcout de calcul devient 
negligeable. Par exemple, cette situation se presente dans un contexte interactif ou la 
solution initiaie <p est prise egale a la solution avant que quelques modifications locales 
5 soient effectuees par I'utilisateur. Dans ce cas, malgre le leger surcout engendre, la 
forme locale de Tequation DSI est preferable car elle est bien plus facile a mettre en 
oeuvre que la forme globale. 



Afin de montrer que I'algorithme propose converge, on peut remarquer que iZ*(y>) 
peut etre exprime sous forme d'une fonction de tp a ayant la forme suivante : 

10 R*{<p) = A-<pl + B-<p a + C 

Les coefficients A, B et C sont independants de tp Q et on a : 

r a = £ K*)o> a (*)) 2 +:E>f 

k€N(a) i 

5 = 2 - { £ W*) *"(*)■ £ «"(*)-^) + £{rf -Q?>} 
y eta 

Le minimum du critere global Bffa) est atteint pour la valeur <p a = -B/{2A) qui 
1 5 correspond precisement a la valeur fournie par i'equation DSI(a). Ceci nous permet de 
conciure que notre algorithme converge car, a chaque iteration du processus iteratif, la 
valeur de la fonction positive ou nulie R?(<p) decroit : 

2«5 Prise en compte de donnees imprecises 

Dans ce qui suit, nous introduisons deux types de contraintes lineaires qui sont d'une 
20 extreme importance en modelisation geometrique ainsi qu'on le verra au paragraphe 
2.7.6. L'objectif de cette section est de prendre en compte des donnees imprecises 
concernant les valeurs inconnues {tp\ : A 6 I}. De telles donnees sont supposees 
associees a des coefficients de certitude w 2 6 R + qui sont en fait des coefficients de 
ponderation positife. 

25 2.5.1 Cas d'une donnee ponctuelle imprecise 

Soit A 6 / un indice de noeud pour lequel la valeur inconnue <p x est supposee proche 
d'une donnee incertaine £ A avec un coefficient de certitude donne egal a vp* : 
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ou correspond a la condition suivante : <p x : 

™l • I^a - <P\\ 2 minimum 

Une telle condition est aisement prise en compte par la methode DSI sous forme d'une 
i eme contrainte avec : 



Akh 



_ f +1 si k = h = A , _ f ¥> x si Jfe = A 
\ 0 sinon 1 | 0 sinon 

et £>** = j 



^ si it =s /i = A 
0 sinon 

Les coefficients Tf , Qf et 7? correspondants utilises dans liquation locale DSI(a) sont 
1 0 alors definis par : 



rf = 0 Qf = J «a b *a « a 
[ 0 sino 

et 7f = | 



= A 
smon 



G7 A si a = A 
0 sinon 

1 5 2.5.2 Cas d'une donnee differentielle imprecise 

Soit A 6 J un index de noeud pour lequel la valeut inconnue <p x est supposee liee a une 
autre valeur <p„ correspondant a un autre index p £ A. Supposons que ce lien ait la 
forme suivante ou A A<1 est une valeur donnee : 

20 (*>#.- Pa) * 

Cette relation est supposee imprecise et est interpretee de la facon suivante ou w? est 
un coefficient de certitude : " 

. ^ • \(v>l - <p\) - A v | 2 minimum 

Une telle condition peut etre facilement prise en compte par la methode DSI en tant 
25 que i eme contrainte avec : 



' 4 A = { 

30 



f -1 si k = h = A 

-1 sifc = /»= M r + a Am si Jfe = A 

+1 si* = Aet/» = M B? = l - A A „ sifc = M 

+1 si& = ^etA = A ( o sinon 

0 sinon 
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r 2 w lti si fc = A = A 



et Df h = < 



2 w a m sik = h=n 



sinon 



Les coefficients Tf , Qf et 7? correspondants utilises dans l'equation locale DSI(a) sont 
alors definis par : 



I? = 



-s?V ' si a = A 
-^•<Px sia = n 
0 sinon 



Q? = 



A A/i si a = A 
sinon 



t*J M si a = A 
0 sinon 

2.5.3 Choix des coefficients de certitude 

Considerons le terme M* de l'equation DSI(a) : 

* Ml = M a 4-Ti +7? + ••■ 
15 avec:M a = £ M (fe) (v a (A)) 2 

Dans les deux exemples presentes ci-dessus, les termes jf sont soit egaux a zero, soit 
correspondent au coefficient de certitude : 

7? = *? 

Cette remarque nous suggere de choisir pour wj un pourcentage donne p { de M a : 
20 &i = Pi • M a Pi >Q 

Le terme M* de liquation DSI(a) est alors egal a : 

Mi = M a .(l+ Pl (a)+p 2 (a) + ...) 
ou p { {<x) est soit egal a un pourcentage donne p f ou egal a zero, 

2.6 Choix des coefficients de ponderation 

25 Le choix des coefficients de ponderation {v a (k)} et <>(*)} est completenient libre ex- 
cepte que les coefficients doivent etre positife ou egaux a zero. Dans ce qui suit 

nous donnons un exemple de choix de ces deux families de coefficients. 
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2.6.1 Choix de coefficients {v a (k)} : 
ponderation harmonique 

Soit A(i) T'orbite" de k definie de la facon suivante : 

A(*) = #(*)-{*} 

Soit |A(*)| le nombre d'elements de A(fc). Nous dirons que nous utilisons une ponderation 
harmontqut si les coefficients {v a {k)} sont choisis de la facon suivante : 

\ 1 si a € A(fc) 
Les fonctions harmoniques ont la propriete caracteristique d'etre egales en tout point 

or iM P I OPre moyenne sur un cercle centr « en ce Point- Comme on peut le voir 
0 si (pfc est egal a la moyenne des valeurs de <p a entouraot le noeud Jb' 
c est pourquoi, par analogie avec les fonctions harmoniques, nous proposons d'appeller 
harmoniques les coefficients {t/ a (fc)} presentes dans cette section. 



Si nous developpons l'equation DSI(a) avec ces coefficients, nous obtenons la for- 
1 5 mule suivante qui peut etre facilement transcrite dans un langage de programmation : 



20 



25 



DSI(a) : 



<Pa = 



jjj-f M(a)-|A(c)[.( £ V*) 

- £ "(*)•{(£ Vfi)-UHk)\. Vi } 

*6A(a) ''WW ' > 



avec 



K = p(aHA(c)| 2 + £ M + £ T ? 

*GA(a) i 



2.6.2 Choix des coefficients {ft(k)} 

Les coefficients {^k)} ont ete introduite afin de moduler Iocalement le lissage de 
1 interpolation. En l'absence de raison particuliere, nous suggerons d'utiliser une 
ponderation uniforme : 

fi(k) = i vfcen 

Parmi toutes les ponderations non uniformes, la ponderation suivante semble parti- 
donn7e? ent mt * teSSaDte S ' U est necessaire d'avoir une fonction lisse au voisinage des 



30 



M(*) 



_ f 1 si k€l 
\ m > 1 si keL 
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2.7 Application a la 

Modelisation Geometrique 

Dans cette section, nous nous proposons de montrer comment utiliser la raethode DSI 
afin d'estimer la surface triangulee 5 elle meme pour laquelle on suppose disposer des 
5 informations partielles suivantes : 

• la position exacte de quelques sommets de triangles, 

• la position approchee de quelques sommets de triangles, 

• des contraintes vectorielles precisant la forme de S, 

• etc . . . 

1 0 Pour cela, on est amene a definir $ a comme le vecteur courant joignant l'origine de 
R au noeud courant de Q et nous noterons <p x , <pV et <p* ses composantes sur la base 
orthonormee {ok, oy, 6z} de R 3 : 

<P = {<P*,<P y ,<P*) 

2.7.1 Definition d'un critere de rugosite globale R(ip) 

1 5 Suivant notre habitude, l'ensemble Q des noeuds de la grille est partage en deux sous- 
ensembles I et L 



<Pi est inconnu V » € / 
fit est connu V I EL 



Les noeuds / 6 L dont la position 0, 6 R 3 est connue seront appeles noeud de 
20 controle et l'on se propose de determiner la position des points restants i € / de telle 
facon que le critere local 2?(£|&) defini ci-dessous et base sur un ensemble donne de 
ponderateurs {v a (k)} soit aussi petit que possible : 

Rim =\\ e «*(*)' vu 2 

<*eiv(*) 

Ce critere constitue la forme vectorielle du critere DSI local. Le critere DSI global 
25 associe R&) est alors defini de la facon suivante ou les {/*(*)} sont des coefficients de 
ponderation non negatifs donnes : 
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Si Ton definit R(<p x ) , R{f^) et R{<f?) ainsi qu'il a ete fait dans les paragraphes precedents : 

k a6iV(fc) 

RW) = E"(t)-W) avec : R(<p"\k) = | X) »*(*) ' fl ? 

oeiV(Jt) 

R{f z ) = !>(*)• avec:i2(^|fc) = | £ • ^ | 2 

* a€JV(fc) 

5 alors, en utilisant le theorems de Pythagore, on verifie facilement que Ton a : 

Rffl = R{<f) + R{<P y ) + R(<P z ) 

2.7.2 Prise en compte de contraintes lineaires au sens des moindres 
carres 

Considerons les matrices suivantes dimensionnees de telle fagon que iV represente le 
10 nombre total de noeuds de Tensemble Q : 

est uae naatrice diagonale non negative de taille N x N 

V v = x, y, z : < [AfJ est une matrice carree de taille N x N 

h est une matrice colonne de taille N 

Si nous voulons que les composantes , et <p z satisfassent les contraintes suivantes 

alors on est araene a introduire le critere suivant 

p(^) = p(^) + p(^) + p(^) 

p(^) = £ IIM-J-^-^llir 
rf*") = £ II -a* III* 



20 



avec : < 



2.7.3 Solution du probleme 

Parmi toutes les surfaces $ interpolant les donnees {$(1) = (p x : I 6 I}, nous nous 
25 proposons de selectionner celle qui minimise le critere suivant : 

R*m = Rtjp) + p(0) 



2652180 



45 



10 



II est facile de verifier que l'on peut toujours decomposer cette formule de la facpn 

SU.lV3»ntG * 

avec: < = W) + p(^) 

Compte tenu de independence des composantes yr*, <p y et on en conclut que : 



R*(<p) minimum 



iZ*(y> r ) minimum 
i?*(y? y ) minimum 
if*(y>*) minimum 



La minimisation du critere vectoriel DSI est done equivalente a la minimisation des 
trois criteres DSI correspondant aux trois composantes du vecteur courant 



15 



2.7.4 Remarque 

Si la surface S est composee de plusieurs morceaux disjoints, alors le theoreme presente 
au paragraphe 2.3 nous garantit que notre probleme admet une solution unique des 
lors qu'il y a au moins un sommet de facette ayant une position connue sur chaque 
morceau. 



2.7.5 Utilisation interactive 

En modelisation geometrique, 1'objectif est de definir interactivement la position dans 
20 Pespace du graphe G associe a Pensemble des noeuds fl. Pour ce faire, on suppose 
connue une forme initiale et Putilisateur deplace quelques noeuds et modifie interac- 
tivement un certain nombre de contraintes (noeuds de controles et donnees imprecises); 
ces modifications sont alors propagees aux noeuds i € I" en utilisant la forme vectorielle 
de la methode DSL 

25 Nous avons montre dans les paragraphes precedents que ce probleme se decomposait 

en trois sous-problemes correspondant aux trois coordonnees (x,y, z) et la meilleure 
fagon de resoudre ces trois problemes est d'utiliser Palgorithme iteratif associe a la 
forme locale des equations DSI correspondantes. 



D'une etape a Pautre du processus de modelisation on peut modifier non seulement 
30 la position des coordonnees des noeuds de contrdle mais Pensemble L lui mime peut 
etre change si necessaire; par exemple, si nous voulons appliquer DSI seulement a un 
sous ensemble T de noeuds, alors nous devons choisir L de telle fa$on que L D f ou f 
designe Pensemble coraplementaire de T dans Q. 

De plus, compte tenu de la decomposition de notre probleme en trois sous-problemes 
3 5 resolus independemment et correspondant aux trois coordonnees, on peut, si necessaire, 
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definir des sous ensembles differents L r , L v et L z pour chacune de ces coordonnees; 
par exempie, si les coordonnees (x,y) ne doivent pas etre modifiees, alors il suffit de 
prendre L x et L y egaux a 1'ensemble Q tout entier. 

2.7.6 Prise en compte de donnees geometriques imprecises 

5 Initialement, la methode DSI a ete congue pour modeliser les surfaces complexes ren- 
contrees dans le domaine des sciences naturelles comme par exempie la geologie et la 
biologie. Dans ce cas, on doit generalement prendre en compte de nombreuses donnees 
impresices comme par exempie : 

f • des noeuds de controle flous 
10 I • des donnees vectorielles floues 

Noeuds de controle flous 

Par definition, nous appellerons "noeud de controle flou" tout noeud A 6 / dont la 
position £ A est connue comme etant approximativement egale a un vecteur donne <p x 
et ceci avec un certain degre de certitude w\. En se referant au paragraphe 2.5.1, on 
1 5 voit que Ton peut prendre en compte une telle information imprecise a Taide de la 
methode DSI en introduisant les contraintes suivantes dans les criteres #*(^ z ), iT r (v? y ) 
et R*{<p z ) : 

( <P\ a 9 A 

20 m^^A 
Contraintes vectorielles floues 

Dans de nombreuses applications, ii est necessaire de controler ie vecteur (£ M - <p x ) 
joignant deux noeuds £ A et £ M d'une surface. Par definition, une telle contrainte sera 
appellee "contrainte vectorielle floue" si elle doit etre satisfaite avec un certain degre 
25 de certitude Quelques exeraples de contraintes de ce type sont presentes sur la 

Figure 2. 

En pratique, au moins un des deux points <p x ou a une position inconnue, par 
exempie celui correspondant a A 6 I, et parmi les cas importants, on peut citer les 
deux suivants : 

30 • Le premier cas correspond aux Figures 2 .a et 2 ,b ou Ton desire controler la 
forme d'un objet de telle fagon que {Jpp — <p\) soit egai a un vecteur donne & Xfi . 

• Le second cas correspond a la Figure 2 .c ou Ton desire controler la forme d'un 
objet de telle fagon que chaque - <p x ) soit orthogonal a un vecteur donne V. 
Cette situation se presente lorsque i'on veut controler le vecteur normal a une 
35 surface; dans ce cas, si V est le vecteur normal au noeud A de la surface, il est 
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judicieux de controler son orthogonalite avec chacun de ces vecteurs (5 -5.) 
pour tout fi € N(X). " 

Si Ton se reporte au paragraphe 2.5.2, le premier cas est facilement pris en compte 
par la methode DSI en introduisant les contraintes suivantes dans les criteres iF^y?*), 
5 R?(<p y ) et avec un coefficient de ponderation egal a wj^ : 

Pour le second cas, le probleme n'est pas tres different car on a : 
10 -<Px)-V~Q 

- <pf) ~ -{vy ■ (y» - + <p' x )}/v* 
Ctp£ - vtt - -{V - - ft) + vy ■ fry - <p»)}/v> 

Dans cecas, nous suggerons de prendre les valeurs des termes(vc^,^,^)et (y? x , <p y , <p z ) 
1 5 intervenant dans les meinbres de droite de ces equations egales aux valeurs correspon- 
dantes obtenues a l'iteration precedente; pour se ramener au cas precedent, il suffit 
alors de poser : 

20 A *m = -{V '-(!£-!P5) + V» -M-rt)}/V 

Remarque 

Dans la prise en compte des donnees imprecises presentee ci-dessus, nous n'avons utilise 
qu'un seul coefficient de certitude wj^ commun aux trois composantes de ces con- 
traintes. En pratique, ceci n'est pas une obligation et, compte tenu de I'independance 
25 des trois composantes, il est toujours possible d'utiliser trois coefficients de certitude 
differents erg, et wf^ correspondant a ces trois composantes. 

2,8 Generalisations 

Avant de conclure, mentionnons quelques generalisations evidentes : 

• La methode DSI peut etre utilisee avec n'importe quels types de facettes polyg- 
30 onales. En pratique, nous conseillons d'utiliser des triangles car ce sont les poly- 
gons les plus simples et les plus faciles a manipuler dans un ordinateur. 

• La methode DSI peut etre utilisee pour estimer des courbes polygonales dans un 
espace 3D; dans ce cas, les triangles sont remplaces par des segments. 
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2.9 CONCLUSION 

La forme locale de ['equation DSI constitue un outil simple et puissant pour modeliser 
les surfaces complexes que i'on peut rencontrer par exemple en geologie ou en biologic 
En particulier : 

• Le maillage utilise pour modeliser une surface peut etre completement libre et 
il est possible d'utiliser des algorithmes automatiques pour les construire; par 
exemple, si le maillage est compose de triangles, on peut utiliser des algorithmes 
derives de la methode de Delaunay (voir ref. [2] et [5]). De plus, la taille du 
maillage peut facilement etre ajustee iocalement en fonction de la complexity de 
la surface, 

• La methode DSI permet de prendre en compte un tres grand nombre de donnees 
precises ou imprecises. Cette possibility est particulierement interessante dans 
le domaine des sciences naturelles ou Tobjectif n'est pas de generer des surfaces 
esthetiques mais d'ajuster ces surfaces aux donnees precises et imprecises de fagon 
coherente et efficace. 

• L'algorithme derive de Tequation DSI locale est suffisamment rapide pour per- 
mettre une utilisation interactive. 

Nous n'avons rien dit sur la representation' des facettes elles-memes. Selon les appli- 
cations, on peut prendre des facettes planes mais, si cela est necessaire, il est toujours 
possible de representer ces facettes par des morceaux de surfaces non planes interpolant 
les noeuds du maillage (voir ref. [4]). 
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RE VEND IC AT IONS 

1. Proc6d6 d'obtention d'une mod61isation d'une 
surface (S) representative notamment de 1' interface entre 
deux zones de natures ou de propri6t6s diff6rentes dans 
un corps & trois dimensions tel qu'une formation 
g6ologique ou un corps vivant, caract6ris§ en ce qu'il 
comprend les 6tapes consistant hi 

- obtenir au moyen d'appareils de mesure un 
ensemble de donn6es g6om6triques relatives k la surface 
et associ6es respectivement 4 des points donn6s de ladite 
surface; 

- effectuer un maillage de la surface de telle 

sorte que 1' ensemble desdits points donn6s soit un sous- 

— ^ 

ensemble des noeuds (C^j) du maillage; 

- m6moriser h une adresse donn6e d'une m6moire, 
sp6cifique k chaque noeud ( Cffc ) du maillage, les 
informations suivantes : 

x ^ l 

* les coordonn&es (Cp k , C^) du noeud consid6r6, 

* le nombre (nb-sat) de noeuds satellites du 
noeud consider 6 , 

* une information (sat) permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives h ces derniers, 

* le cas 6ch£ant des donn6es g6om6triques 
(const) associ6es audit noeud consid6r§, 

- pour chaque noeud du maillage f dfefinir un indice 
de rugosity local (R(C^|k)) obtenu A partir d'une somme 
pond6r6e des coordonn6es actuelles du noeud et de ses 
satellites, 

- d6finir la somme (R* (Op) ) d'un indice de rugosit6 
global (R«f)) f constitue par l'addition des indices de 
rugosit§ locaux associ6s & chaque noeud f et d'un indice 
global (^>(^)) de violation desdites donn6es gfeom6triques ; 

- ajuster de fagon iterative les coordonn§es de 
chacun des noeuds dont on ne connalt pas les coordonn6es 
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pr6cises, en utilisant k chaque fois 1' addition d'une 
combinaison pond6r6e des coordonn6es actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donnfees g6om6triques associ6es 
audit noeud f de mani^re k amener ladite somme k un 
minimum; 



- cr6er une representation de la surface k partir 



2. Proc6d6 selon la revendication l f caract6ris6 
en ce qu'on effectue un maillage triangulaire de la 
surface, 

3. Proc6d6 selon l'une des revendications 1 et 2, 
caract6ris6 en ce que l'etape d' ajustement comprend trois 
sous-etapes distinctes mettant en jeu respectivement les 
trois coordonn§es spatiales (^*, ^tflc) des noeuds et des 
satellites . 

4. Procede selon l'une des revendications 1 k 3, 

caract6rise en ce que les donn6es g§om6triques 

— =*> 

comprennent des donn6es de coordonn6es (<|>k ) des noeuds . 

5. Proc6d6 selon l'une des revendications 1 k 4, 
caract6ris6 en ce que les donn6es g6om6triques 
comprennent des donn6es de vecteurs (C^-tfk) entre deux 
noeuds donn6s . 

6. Proced6 selon l'une des revendications 1 k 5, 
caract6ris6 en ce que les donn6es g*6om6triques 
comprennent des donn6es de vecteur normal (V) k la 
surface k mod61iser. 

7. Proc6d6 selon l'une des revendications 4 & 6 r 
caract6ris£ en ce qu'au moins certaines desdites donn6es 
g6om6triques sont assorties d'un coefficient (TO^r^ty* 
repr6sentatif du degr6 de certitude avec lequel ces 
donn6es sont connues et en ce que l'indice de violation 
des donn6es int6gre lesdits coefficients de certitude. 

8. Proc6d§ selon l'une des revendications 
pr6c6dentes, caract6ris6 en ce qu'il comprend en outre 
les Stapes consistant k mesurer des propri6t6s du corps k 



des coordonn§es ajust§es 




chacun des noeuds. 
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trois dimensions dans la region d'au moins certains 
noeuds du maillage et k m6moriser les mesures obtenues 
(info) aux adresses sp£cifiques aux noeuds consid6res. 

9. Proc6d6 selon l'une des revendications 
5 pr6c6dentes f caract6ris6 en ce que r dans chaque 

coefficient de rugosite local , le rapport entre le poids 
(v° (k) ; a=k) du noeud concern^ et le poids (v a (k) ; ae 
A>{tf) d'un satellite quelconque dudit noeud est 6gal au 
nombre (nk ) de satellites dudit noeud, affects d*un signe 
negatif. 

10. Precede selon la revendication 9, caract6rise 
en ce que ladite combinaison pond6r6e utilise les m§mes 
poids . 

11. Procede selon l'une quelconque des 
15 revendications pr§cedentes f caract6ris6 en ce que la 

somme des coefficients de rugosite locaux est une somme 
pond6ree. 

12. Procede selon la revendication 11, rattach6e k 
la revendication 10 , caracterise en ce que ladite 

20 combinaison pond6r6e utilise 6galement les poids (y(k)) 
associes aux coefficients de rugosite locaux. 

13. Procede selon l'une des revendications 
prec6dentes, caracterise en ce que la surface est une 
surface de separation entre deux milieux de propriet6s ou 

25 de natures diff6rentes dans une formation g£ologique et 
en ce qu'on utilise la representation de la surface pour 
optimiser 1 1 exploration de richesses souterraines de 
ladite formation. 

14. Proc6de selon l'une des revendications 1 k 12, 
caracterise en ce que la surface est representative des 

30 

variations d'un propriety d'un corps k deux dimensions au 
voisinage d'un plan traversant ledit corps, ladite 
propriety 6tant representee par la dimension restante, et 
en ce qu'on utilise la representation de la surface pour 
optimiser 1 1 exploration ou 1 1 exploitation de richesses 

35 

souterraines de ladite formation. 
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15. Procede selon l'une des revendications 1 k 14, 
caracteris6 en ce que la representation de la surface est 
une representation graphique sur un support plan. 

16. Proc6de selon l'une des revendications 1 k 14, 
caract6ris6 en ce que la representation de la surface est 
une maquette en trois dimensions. 

17. Dispositif pour l'obtention d'une modeiisation 
d'une surface (S) representative notamment de l 1 interface 
entre deux zones de natures ou de propri6tes diff6rentes 
dans un corps k trois dimensions tel qu'une formation 
geologique ou un corps d'etre vivant, caracterise en ce 
qu'il comprend : 

- des moyens de mesure un ensemble de donnees 
geometriques relatives k la surface et associ6es 
respectivement k des points donn6s de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donnes soit un sous-ensemble des noeuds du maillage; 

- une memoire contenant k des adresses donnees, 
sp6cifiques cl chaque noeud du maillage, les informations 
suivantes : 



* le nombre (nb-sat) de noeuds satellites du 
noeud considere, 

* une information (sat) permettant d'acc6der aux 
adresses specifiques desdits noeuds satellites et par 
suite aux informations relatives k ces deirniers, 

* le cas ech6ant des donnees geometriques 
(const) associees audit noeud consid6r6, 

- des moyens de calcul pour ajuster de fa<pon 
iterative les coordonnees de chacun des noeuds dont on ne 
connait pas les coordonnees precises, en utilisant k 
chaque fois l 1 addition d'une combinaison pond6r6e des 
coordonnees actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
donnees g6om6triques associ6es audit noeud, de mani^re k 



* les coordonnees 
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amener k un minimum une somme (R*(C|?)) d'un indice de 
rugosit6 global (R(Cf)), constitu6 par 1* addition 
d 1 indices de rugosit6 locaux (R( Cp j k)) respectivement 
associ6s h chacun des noeuds et chacun obtenu & partir 
d'une somme pond6r6e des coordonn6es actuelles dudit 
noeud et de ses satellites f et d'un indice global (^(^)) 
de violation desdites donn6es g£om6triques; et 

- des moyens pour cr6er une representation de la 
surface a partir des coordonn6es ajust6es de chacun des 
noeuds . 



2652180 



1/6 




FIG.1 



• noeud de con f role 
onoeud libra 




FIG. 2a 





'CABINET REGIMBEAU 
MARTIN, SCHRKVIPE 
WARCOM, AHNER 

3UPP ?lk ~ 

certify conforme I I'origuial 



FIG. 2c 




FIG. 2d 



2/ 6 



2652180 




FIG.3 




«rtffie* coaform3 a 1'original 



2652180 



3/6 




Cabinet regimbeau 
MARTIN, SCHRMPF 
WARG01N, AHNifi 

certify catena I 




2652180 



4/6 



sal" 
(=L) 





adresse(t) 




tableau 


adresse(t+f 


A*(j 2 )- 


de . 






n peases 








adresseCt+n^ 





A(j 2 ) 



AOnk) 



FIG_6 



sat- 
C«il) 



adresseCL 2 ) ^Jg^ 



A*G2) 



adresse(ink) 



next 
code fin) 



A*(Jnk)- 



AGO 

AG 2 ) 



AOnk) 



CABINET REGIMBEAU 
MARTIN, SCHRfMPF 
WARCOiN, AKNER 

OUPLeGATA 

csrtifie conforms & I'origiaal 



FIG.7 



2652180 



5/6 




FIG.8 



Afome A(k) 



noyau 


nb-sat 


sat" 


consf 


info 






(-adresseQ] 





Cabinet regimbeau 
MARTIN. SCHRMPF 
WARCOM, AKNER 

certify canfanns a I'original 



Uste des conlraintes 



adresse (Ci) 



adresse(C2) 



(*adr 


esseC2) 


contrainte 


i 




nex^ 
(= adresseC3) 


conhrainte 


i 


* 

- 


next" 
(» code fin) 


contrainte 




FIG. 9 



6/6 



2652180 



A (1) 




AC1) 


(sadressea) 


► 


A (2) 




AC2) 


(*adressea*i) 


p- 








A (N) 




ACN) 


(<*adressea+N-l) 


> 



FIG.10 



SRBINET REGIHBEfttf 
MARTIN, SCHRBJIPF 
WARCOIN, MINER 

certify conforms & I'original 



A*(D 
(»adresseai) 



A*(N) _ 
(=adresseaN) 



A*(2) 




C e aaresse B2) 


^ 



nexh 
(«adressea 2 ) 


Atom 
A(1) 


i 


^ex^ 
(=adresse33) 


Atom 
A(2) I 


i 


next 
(©code fin) 


Atom 
A(N) 



FIG.11 



REPUBLIQUE FRANCAISE 



2652180 



N° d'enregistrement 



INSTITUT NATIONAL 
de la 

PROPRIETE INDUSTRIELLE 



RAPPORT DE RECHERCHE 

etabli sur la base des dernieres revendtcations 
deposees avant le commencement de la recherche 



FR 8912341 
FA 431687 



Categorie 



DOCUMENTS CONSIDERES COMME PERTINENTS 



Citation du document avec indication, en cas de besoin, 
des parties pertmentes 



Revendications 
concern ecs 
de la demande 



COMPUTERS & GRAPHICS, vol. 6, no. 2, 
1982, pages 63-72, Oxford, GB; SINGH, 
SPILLERS: "FEM shape optimization: A 
case for interactive graphics" 
* Page 65, colonne de gauche, ligne 34 
- colonne de droite, ligne 8 * 

EP-A-0 196 651 (HITACHI) 



1-4,8, 
14-17 



DOMAINES TECHNIQUES 
RECHERCHES (Int. OS) 



G 06 F 15/72 



a 
2 



Date (TacMTcmfnt de la recherche 

11-06-1990 



BURGAUD C. 



CATEGORIE DES DOCUMENTS CITES 

X : partlculierement pertinent £ !ui seul 

Y : particulierement pertinent ea combinaisan avec tin 

autre document de ta meme categorie 
A : pertinent a rencontre d'au moins une revendicatlon 

ou arriere-plan technologique general 
O : divulgation non-ecrite 

P : document intercalaire 



T : theorie ou principe a la base de Y invention 

E : document de brevet beneficiant d'une date anterieure 

i la date de depflt et qui n'a ete publie qir'i cette date 

de depot ou qir'i une date posterieuxe. 
D : cite dans la demande 
L : cite pour d'autres ralsons 

& : membre de la meme famille, document correspondant 



